Measurement of the mass difference between t and t quarks 
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We present a direct measurement of the mass difference between t and t quarks using tt candidate 
events in the lepton+jets channel, collected with the CDF II detector at Fermilab's 1.96 TeV Teva- 
tron pp CoUider. We make an event by event estimate of the mass difference to construct templates 
for top quark pair signal events and background events. The resulting mass difference distribu- 
tion of data is compared to templates of signals and background using a maximum likelihood fit. 
From a sample corresponding to an integrated luminosity of 5.6 fb~^, we measure a mass difference, 
AMtop ~ Mt — Mt = —3.3 ± 1.4 (stat) ± 1.0 (syst) GeV/c^, approximately two standard deviations 
away from the CPT hypothesis of zero mass difference. 

PACS numbers: 14.65. Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 
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Discrete symmetries reflecting the invariance under 
discrete transformations, such as charge conjugation (C), 
space reflection or parity (P), and time reversal (T), are 
not always exact. Examples include the C and P sym- 
metries and their CP combination, which are violated 
by the weak interactions CPT symmetry, which re- 
flects the invariance under the combined operation of C, 
P, and T transformations, has not been found to be vi- 
olated in any experiment so far 0, However, it is 
important to examine the possibility of CPT violation 
in all sectors of the standard model (SM), as there are 
well-motivated extensions of the SM allowing for CPT 
symmetry breaking In the CPT theorem, particle 
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and anti-particle masses must be identical; thus, a mass 
difference between particle and its anti-particle would in- 
dicate a violation of CPT. The mass equality has been 
verified to high precision for leptons and hadrons, but 
not for quarks. With the exception of the top quark, it 
is impossible to measure quark masses directly, because 
a newly created quark dresses itself with other quarks 
and gluons to form a hadron, and hadron masses yield, 
at best, only rough estimates of the quark mass. The 
top quark is by far the most massive quark and, with 
lifetime of the order of 10"^'' seconds, decays before it 
can hadronize. This allows a precise measurement of the 
mass difference between t and t quarks and provides a 
probe of CPT violation in the quark sector . 

This letter reports a measurement of the mass differ- 
ence (AMtop — Mt — Mf ) between t and t quarks us- 
ing a sample of tt candidates in the lepton-|-jets final 
state. The data correspond to an integrated luminosity 
of 5.6 fb""'^ in proton-antiproton collisions at the Tevatron 
with = 1.96 TeV, collected with the CDF II detec- 
tor ^] . Assuming unitarity of the three-generation CKM 
matrix, t and t quarks decay almost exclusively into a W 
boson and a bottom quark 

The case where one W decays into a charged lepton and 
a neutrino {W^ — )■ iiy or W~ — iD) and the other into 
a pair of jets defines the lepton-t-jets decay channel. The 
electric charge of the lepton (-1 for £ and for i) deter- 
mines the fiavor of top quarks with event reconstruction. 
To select tt candidate events in this channel, we require 
one electron (muon) with Ex > 20 GeV {pr > 20 GeV/c) 
and pseudorapidity \r]\ < 1.1 Q. We also require high 
missing transverse energy Q, > 20 GeV, and at 
least four jets. Jets are reco nstructed with a cone algo- 
rithm 9] with radius R = y/i^v)^ + (Ai?!))^ =0.4. Jets 
originating from b quarks are identified using a secondary 
vertex tagging algorithm [lol |. In order to optimize the 
background reduction process and improve the statistical 
power of the events, we divide the sample of tt candidate 
events into sub-samples with zero, one, and two or more 
6-tagged jets. 

When an event has zero or one 6-tagged jet, we require 
exactly four jets with transverse energy Et > 20 GeV 
and 1 77 1 < 2.0. If an event has two or more 6-jets, three 
jets are required to have Et > 20 GeV and I77I < 2.0, 
and a fourth jet is required to have Et > 12 GeV and 
I77I < 2.4, with no restriction on the total number of 
jets. To reject backgrounds, we require the scalar sum 



of transverse energies in the event, Ht = E, 



lepton 



-^7 



Sfour jets ' to be greater than 250 GeV. 

The primary sources of background events are M^-l-jets 
and QCD niultijet production. Contributions from 
Z-|-jets, diboson, and single top production are expected 
to be small. To estimate the contribution of each process, 
we use a combination of data and Monte Carlo (MC) 



based techniques described in Ref. [11|. For the Z-l-jets, 
diboson, and single top quark events, we normalized MC 
simulation events using their respective theoretical cross 
sections. The QCD multijet background is estimated 



TABLE I: Expected and observed numbers of signal and 
background events assuming ti production cross-section att 
= 7.4 pb and Mtop = 172.5 GeV/c^ 



6-tag 



1 6-tag 



> 2 &-tag 



W+jets 


596 ± 98 


88.3 ± 


23.0 


11.1 ± 


3.6 


QCD muhijet 


95.8 ± 74.4 


14.7 ± 


12.1 


2.4 ± 


3.2 


Z-l-jets 


48.8 ± 9.4 


5.7 ± 


1.3 


0.8 ± 


0.2 


Diboson 


50.1 ± 4.7 


6.6 ± 


0.8 


1.0 ± 


0.2 


Single top 


4.0 ± 0.4 


5.5 ± 


0.5 


2.2 ± 


0.2 


Background 


795± 124 


121 ± 


24 


17.3 ± 


4.8 


tt signal 


426 ± 57 


578 ± 


72 


282 ± 


44 


Expected 


1220 ± 137 


699 ± 


76 


299 ± 


44 


Observed 


1278 


720 


296 



with a data-driven approach. We model I^-l-jets back- 
ground events using MC simulation but the overall rate 
is determined using data after subtracting the rate of all 
the other backgrounds and tt. Table |T] shows the expected 
background composition and the expected number of tt 
events. 

We assume selected events to be tt events in the Icp- 
ton-f-jets channel and reconstruct them to form estima- 
tors of AMtop, using a special purpose kinematic fitter, 
in which we modify the standard fitter [l^ to allow a 
mass difference between t and i. Measured four-vectors 
of jets and lepton are corrected for known effects [l3[, 
and resolutions are assigned. The unclustered transverse 
energy (Ut), which is the sum of all transverse energy in 
the calorimeter that is not associated with the primary 
lepton or one of the leading four jets, is used to calcu- 
late the neutrino transverse momentum. The longitudi- 
nal momentum of the neutrino is a free (unconstrained) 
parameter which is effectively determined by the con- 
straint on the invariant mass of the leptonic W. We then 
define a kinematic fit having a free parameter drrireco , 



'^i=iAjetsij>T 
-'k—x,y 



ijit 



Pi ^ 



y, (TT^'^'^ — TTmeas\2 / 2 



(Mj, - Mwf/Tl^ + (M,, - Mwf/rl^ 



{Mbn - (Mtop ■ 
{Mm^ - (Mtop 



dnir 
drrir 



,/2)}Vr? 
o/2)}Vr? 



(1) 



where dm™", the drrireco value at the lowest x ^ rep- 
resents the reconstructed mass difference between the 
hadronic and leptonic top decay (Mf,jj — M^j,). In this 
formulation, the first term constrains the pT of the lep- 
ton and four leading jets to their measured values within 
their uncertainties (ui); the second term does the same 
for both transverse components x and y of the unclus- 
tered transverse energy. In the remaining four terms, 
the quantities Mj ^ , M^^ , Mf,^ j , and yibiv refer to the in- 
variant masses of the four vector sum of the particles 
denoted in the subscripts. M-w and Mtop are the masses 
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of the W boson (80.4 GeV/c^) Q and the average of t 
and t quark masses (172.5 GeV/c^V close to the current 
best experimental determination [l3|, respectively. Tw 
(2.1 GeV/c2) and Ft (1.5 GeV/c^) are the total widths 
of the W boson and the t quark l| . We assume that the 
total widths of the t and t quarks are equal. Determining 
the reconstructed mass difference of t and f, Amioco, re- 
quires the identification of the flavor (t versus t) , and this 
is done using the electric charge of the lepton {Qiepton) , 

defining Amrcco = -Qiepton X dm™". 

The use of different detector components and the dif- 
ferent resolutions of the measured values for jet, lepton, 
and unclustered energy, make the reconstructed mass dis- 
tribution of hadronic top quarks differ from that of lep- 
tonic top quarks. Because the sign of Amj-cco depends 
on the lepton charge, Amj-cco distributions for the posi- 
tive and negative lepton events are different. We divide 
the sample into six sub-samples, two samples with posi- 
tively and negatively charged leptons for each of 6-tag, 
1 6-tag, and 2 6-tag samples. 

With the assumption that the leading four jets in the 
event come from the four final quarks at the hard scat- 
tering level, there are 12, 6, and 2 possible assignments 
of jets to quarks for 6-tag, 1 6-tag, and 2 &-tag re- 
spectively. The minimization of is performed for each 
jet-to-parton assignment, and Amrcco is taken from the 
assignment that yields the lowest iXmin)- Events with 
Xmin > 9.0 {Xmin > ^.0) are removed from the sample 
to reject poorly reconstructed events for 6-tagged (zero 
&-tagged) events. To increase the statistical power of 
the measurement, we employ an additional observable 

(2) w 

Amreco from the assignment that yields the 2 lowest 
X^- Although it has a poorer sensitivity, Amrcco pro- 
vides additional information on AMtop and improves the 
statistical uncertainty by approximately 10 %. 

Using MADGRAPH [l^, we generate ti signal samples 
with AMtop between -20 GeV/c^ and 20 GeV/c^ us- 
ing almost 2 GeV/c^ step size, where we take the aver- 
age mass value of t and t to be Mtop = 172.5 GeV/c^. 
Parton showering of the signal events is simulated with 
PYTHIA [l^, and the CDF detector is simulated using a 
GEANT-based software package [itI] . 

We estimate the probability density functions (p.d.f.s) 
of signal and background templates using the kernel den- 
sity estimation (KDE) [H, [l3| . For the AMtop measure- 

(2) 

ment with two observables (Amreco and Amreco), we use 
the two dimensional KDE that accounts for the correla- 
tion between them. First, at discrete values of AMtop 
from —20 GeV/c^ to 20 GeV/c^, we estimate the p.d.f.s 
for the observables from above-mentioned tt MC sam- 
ples. We interpolate the MC distributions to find p.d.f.s 
for arbitrary values of AMtop using the local polynomial 
smoothing method Q . We fit the signal and background 
p.d.f.s to the measured distributions of the observables in 
the data using an unbinncd maximum likelihood fit [2ll | , 
where we minimize the negative logarithm of the likeli- 
hood with MINUIT tea. Likelihoods are built for each of 



TABLE II: Summary of systematic uncertainties on AMtop. 



Source 


Uncertainty (GeV/c'^) 




0.7 




0.4 




0.2 


Pctrtoii distribution functions 


0.1 


h-\ft pTipfcfV apfllp 


0.1 


Background shape 


0.2 


Gluon fusion fraction 


0.1 


Initial and final state radiation 


0.1 


Monte Carlo statistics 


0.1 


Lepton energy scale 


0.1 


Multiple hadron interaction 


0.4 


Color reconnection 


0.2 


Total systematic uncertainty 


1.0 



six sub-samples separately, and an overall likelihood is 
then obtained by multiplying them together. We eval- 
uate the statistical uncertainty on AMtop by searching 
for the points where the negative logarithm of the likeli- 
hood exceeds the minimum by 0.5. Refs. [3 113 provide 
detailed information about this technique. 

We test the fitting procedure using 3000 MC pseudo 
experiments (PEs) for each of 11 equally spaced AMtop 
values ranging from —10 GeV/c^ to 10 GeV/c^. The dis- 
tributions of the average residual of measured AMtop (de- 
viation from the input AMtop) for simulated experi- 
ments is consistent with zero. However the width of the 
pull (the ratio of the residual to the uncertainty reported 
by minuit) is 4 % greater than unity. We therefore in- 
crease the measured uncertainty by 4 %. 

We examine a variety of systematic effects that could 
change the measurement by comparing results from PEs 
in which we vary relevant systematic parameters within 
their uncertainties. All systematic uncertainties are sum- 
marized in Table |lll The dominant source of systematic 
uncertainty is the signal modeling, which we estimate 
using PEs with events generated with MADGRAPH and 
PYTHIA. We also estimate a parton showering uncer- 
tainty by applying different showering models (PYTHiA 
and HERWIG (2JI) to a sample generated with ALP- 
GEN [2^ . We address a possible difference in the detector 
response between b and b jets by comparing data and MC 
simulation events [2^. We add a systematic uncertainty 
due to multiple hadron interactions to account for the 
fact that the average number of interactions in our MC 
samples is not exactly equal to the number observed in 
the data. The jet energy scale (JES), the dominant un- 
certainty in most of the top quark mass measurements, is 
partially canceled in the measurement of the mass differ- 
ence. Therefore JES contributes only a small uncertainty 
to this measurement. Other sources of systematic effects, 
including uncertainties in parton distribution functions, 
gluon radiation, background shape and normalization, 
lepton energy scale, and color reconnection [1^, H^l, give 
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small contributions. The total systematic uncertainty of 
1.0 GeV/c^ is derived from a quadrature sum of the listed 
uncertainties. 

The likelihood fit to the data returns a mass difference 

AMtop = -3.3 ± 1.4 (stat) ± 1.0 (syst) GeV/c^ 
= -3.3 ± 1.7 GeV/c2. 

Figure [T] shows the measured distributions of the ob- 
servables used for the AMtop measurement overlaid with 
density estimates using tt signal events with AMtop 
— ~4 GeV/c^ and GeV/c^ and the full background 
model. The choice of AMtop = -4 GeV/c^ (sohd 
line) gives better agreement with the data than that of 
GeV/c2 (dashed line). 

In conclusion, we examine the mass difference be- 
tween t and t quarks in the lepton+jets channel us- 
ing data corresponding to an integrated luminosity of 
5.6 fb^^ from pp collisions at ^/s — 1.96 TeV. We mea- 
sure the mass difference to be AMtop = Mt — Mj = 
-3.3±1.4(stat)±1.0 (syst) GeV/c^ = -3.3±1.7 GeV/c^. 
This result is consistent with CPT-symmetry expecta- 
tion, AMtop = GeV/c^, with approximately 2a level 
deviations. It is consistent with the recent result from 
the DO Collaboration [2^, but is 2.2 times more precise. 
This is the most precise measurement of the mass differ- 



ence between t and t quarks. 
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